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Abstract 

We study the propagation properties of nonnegative and bounded solutions of the class 
of reaction-diffusion equations with nonlinear fractional diffusion: u t + (— A) s (u m ) = f(u). 
For all < s < 1 and m > m c = (N — 2s) + / 'N , we consider the solution of the initial-value 
problem with initial data having fast decay at infinity and prove that its level sets propagate 
exponentially fast in time, in contradiction to the traveling wave behaviour of the standard 
KPP case, which corresponds to putting s = 1, m = 1 and f(u) = u(l — u). The proof 
of this fact uses as an essential ingredient the recently established decay properties of the 
self-similar solutions of the purely diffusive equation, u t + (— A) s u m = 0. 
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1 Introduction 



We consider the following reaction-diffusion problem 



(1.1) 



u t (x, t) + L s u m {x, t) = f(u) for x E M JV and t > 0, 
u(x, 0) = Uq(x) for x E 



DiV 



where L s = (— A) s is the Fractional Laplacian operator with s E (0, 1). We are interested 
in studying the propagation properties of nonnegative and bounded solutions of this problem 
in the spirit of the Fisher-KPP theory. Therefore, we assume that the reaction term f(u) 
satisfies 

(1.2) / E C^QO, 1]) is a concave function with /(0) = /(l) = 0, f(l) < < /'(0). 

For example we can take f(u) = u(l — u). Our results will depend on the parameters m and 
s, according to the ranges m c < m < mi, mi < m < 1, and m > 1, where 

(JV-2a)+ iV 



m c = — , mi 



N 1 N + 2s 

1.1. Perspective. The traveling wave behavior. The problem with standard diffusion 
goes back to the work of Kolmogorov, Petrovskii and Piskunov, see [20], that presents the 
most simple reaction-diffusion equation concerning the concentration u of a single substance 
in one spatial dimension, 

(1.3) d t u = Du xx + f(u). 

The choice f(u) = u(l — u) yields Fisher's equation [18] that was originally used to describe 
the spreading of biological populations. The celebrated result says that the long-time behavior 
of any solution of (1.3), with suitable data < uo(x) < 1 that decay fast at infinity, resembles 
a traveling wave with a definite speed. When considering equation (1.3) in dimensions N > 1, 
the problem becomes 



(1.4) u t -Au = f(u) in(0,+oo)x 



which corresponds to (1.1) in the case when L s = —A, the standard Laplacian. This case has 
been studied by Aronson and Weinberger in [31 H], where they prove the following result. 

Theorem AW. Let u be a solution of (1.4) with uq ^ compactly supported in ~K N and 
satisfying < n (-) < 1. Let c* = 2 A // / (0). Then, 

1. if c > c*, then u(x, t) — > uniformly in {\x\ > ct} as t — > oo. 

2. if c < c*, then u(x, t) — > 1 uniformly in {\x\ < ct} as t oo. 



In addition, problem (1.4) admits planar traveling wave solutions connecting and 1, that is, 
solutions of the form u(x, t) = <p(x ■ e + ct) with 

-(f)" + C (f>' = /((/>) in M, 0(-oo) = 0, 0(+oo) = 1. 
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This asymptotic traveling-wave behavior has been generalized in many interesting ways. Of 
concern here is the consideration of nonlinear diffusion. De Pablo and Vazquez study in 
|16j the existence of traveling wave solutions and the property of finite propagation for the 
reaction-diffusion equation 

ut = (u m ) xx + Xu n (l - u), (x, t)eix (0, oo) 

with m > 1, A>0, n G M and u = u(x,t) > 0. Similar results hold also for other slow 
diffusion cases, m > 1, studied by de Pablo and Sanchez (|15j). 

1.2. Non-traveling wave behavior. Departing from these results, King and McCabe 
examined in |19j a case of fast diffusion, namely 

u t = Au m + u(l-u), x€R N ,t>0, 

where (N — 2) + /N < m < 1. They showed that the problem does not admit traveling wave 
solutions. Using a detailed formal analysis, they also showed that level sets of the solutions of 
the initial-value problem with suitable initial data propagate exponentially fast in time. They 
extended the results to all < m < 1. 

On the other hand, and independently, Cabre and Roquejoffre in [l0l IH] studied the case of 
fractional linear diffusion, s G (0, 1) and m = 1, and they concluded in the same vein that there 
is no traveling wave behavior as t — > oo, and indeed the level sets propagate exponentially fast 
in time. This came as a surprise since their problem deals with linear diffusion. 

Motivated by these two examples of break of the asymptotic TW structure, we study here 



the case of a diffusion that is both fractional and nonlinear, namely problem ( 1.1 ) in the range 
s G (0, 1) and m > m c . The initial datum uq(x) : l w — > [0, 1] and satisfies a growth condition 
of the form 



(1.5) < u (x) < C\x\- X( - N ' s ' m \ Vx G 



where the exponent X(N, s, m) is stated explicitly in the different ranges, m c < m < m\ and 
mi < m- In this paper we establish the negative result about traveling wave behaviour, more 
precisely, we prove that an exponential rate of propagation of level sets is true in all cases. 
We also explain the mechanism for it in simple terms: the exponential rate of propagation of 
the level sets of solutions (with initial data having a certain minimum decay for large |x|) is 
a consequence of the power-like decay behaviour of the fundamental solutions of the diffusion 
problem studied in |23| . Therefore, we obtain two main cases in the analysis, m c < m < mi 
and m > mi, depending on that behaviour. 



1.3. Main results. The existence of a unique mild solution of problem (1.1) follows by 
semigroup approach. The mild solution corresponding to an initial datum uq G L}{W n ), 
< uq < 1 is in fact a positive, bounded, strong solution with C 1,a regularity. In the 
Appendix we give a brief discussion of these properties. Let us introduce some notations. 
Once and for all, we put /3 = l/(N(m — 1) + 2s) and 

1 — m ,,, v 1 „,, , 1 + 2(m — l)Bs „,, , 
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The value o\ appears for m c < m < m\ and then o\ > 02- Notice also that 02 < 03 for 
m > 1. Here is the precise statement of our main results for the solutions of the generalized 



KPP problem (1.1) 



Theorem 1.1 Let N > 1, s 6 (0, 1), / satisfying (1.2) and m\ < m < 1. iei u be a solution 



of ( |1.1[ ), where < ito(-) < 1 is measurable, uq 7^ and satisfies 
(1.7) < u (x) < C|x|- (7V+2s) , VxGR^. 

T/ien 

-Z. if cr > 02, then u(x,t) — > uniformly in {\x\ > e' 7 *} as t ^ 00. 
,2. if cr < (72, £/ten u(x,t) — )• 1 uniformly in {\x\ < e' 7 *} as t ^ 00. 



Theorem 1.2 Xei JV > 1, s € (0,1), / satisfying (1.2) and m s < m < m\. Let u be a 



solution of (1.1), where < uo(-) < 1 is measurable, uq 7^ and satisfies 
(1.8) < n (x) < C|a;|- 2s/(1 - m) , VxGl w 

T/ien 

-Z. if cr > o~\, then u(x,t) — > uniformly in {\x\ > e' 7 *} as t ^ 00. 
2. if cr < o\, then u(x,t) — >• 1 uniformly in {\x\ < e°"*} as i — >• 00. 



Theorem 1.3 Xei JV > 1, s G (0,1), / satisfying (1.2) and m > 1. Let u be a solution of 



(1.1), where < uo(-) < 1 is measurable, uq ^ and satisfies 

< u (x) < C\x\-( N+2s \ VxeR N . 

Then 

1. if cr > o"3, i/ien u(x,t) — >• uniformly in {\x\ > e°"*} as t ^ 00. 

2. if cr < <72, i/ien u(x,t) — > 1 uniformly in {\x\ < e at } as t — )■ 00. 



Remarks. In all ranges of parameters m > m c , there appear critical values of <7 with an 
influence on the behavior of the level sets. 

• In the case mi < m < 1, the case cr = 02 is still open. This critical exponent is the same as 
in the case of the linear diffusion m = 1, proved in |llj . 

• In the range m c < m < mi, the case <7 = cii is still open. In particular, for the classical 
case s = 1 and f(u) = n(l — u) we get <7i = which is a critical speed found by King and 
McCabe [19]. In this way, we complete their result with rigorous proofs to all s S (0, 1). 
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Figure 1: Ranges of parameters m and s: we study the cases R\ - 
m < N/(N + 2s)}; R 2 = {s G (0, 1), iV/(iV + 2s) < m < 1}, R 3 



{sG (0,1), (iV-2s)+/iV< 
= {s G (0,1), m > 1} 



• In the case m > 1, we do not cover the entire interval [02,03]. Therefore, we prove that for 
m > 1 the nonlinearity has a different influence on the velocity of propagation. 

• The result of Theorems 1.1 and 1.2 is true also in the case m = mi, where u\ = 02. The 
outline of the proof is the same, but there are a number of additional technical difficulties, 
typical of borderline cases. We have decided to skip the lengthy analysis of this case because 
of the lack of novelty for our intended purpose. 

Our main conclusion is that exponential propagation is shown to be the common occurrence, 
and the existence of traveling wave behavior is reduced to the classical KPP cases mentioned 
at the beginning of this discussion (see dotted line in Figure [T]) . 

As we have already mentioned, one of the motivations of the work was to make clear the 
mechanism that explains the exponential rate of expansion in simple terms, even in this 
situation that is more complicated than |XQ(, lllj. In fact, due to the nonlinearity, the solution 
of the diffusion problems involved in the proofs does not admit an integral representation as the 
case m = 1. Instead, we will use as an essential tool the behavior of the fundamental solution 
of the Fractional Porous Medium Equation, also called Barenblatt solution, recently studied 
in |23j. To be precise, the decay rate of the tail of these solutions as \x\ — > 00 is the essential 
information we use to calculate the rates of expansion. This information is combined with 
more or less usual techniques of linearization and comparison with sub- and super-solutions. 
We also need accurate lower estimates for positive solutions of this latter equation, and a 
further selfsimilar analysis for the linear diffusion problem. 

1.4. Organization of the proofs. In Section |4j under the assumption of initial datum 
with the decay (1.5), we prove convergence to in the outer set {|x| > e at } by constructing a 



super-solution of the linearized problem with reaction term f'(0)u. The arguments hold for a 
larger than the corresponding critical velocity. 



5 



In Section [5] we prove convergence to 1 on the inner sets {|x| < e at } in various steps. We 
only assume < no < 1, uo 7^ 0. We first show that the solution reaches a certain minimum 



profile for positive times, thanks to the analysis of Theorem 1.4 below, we then perform an 
iterative proof the conservation in time of this minimum level, and finally convergence to 1 is 
obtained by constructing a super-solution to the problem satisfied by 1 — u m . Therefore, we 
deal with a problem of the form 

a(x, t) wt(x, t) + L s w(x, t) + bo w(x, t) > 0. 

A suitable choice for constructing the super-solution w is represented by self-similar solutions 
of the form U(x, t) = t a F(\x\t~@ ) of the linear problem 

(1.9) U t + L s U = 

with radial increasing initial data. This motivates us to derive a number of properties of the 



linear diffusion problem (1.9), also known as the Fractional Heat Equation. In particular, we 
need to show that the profile F mentioned above has the same asymptotic behavior as the 
initial data. In order to establish such fact we have to review, Section [6l the properties of the 



fundamental solution of Problem (1.9) 



K s (x,t) = rs/(t-K|s|), f(r) ~ r^ N+2s \ 
We perform a further analysis of the profile / by proving that rf ~ r -(N+2s)_ 

Remark. As a consequence of the exponential propagation of the level sets, we immediately 
obtain the non-existence of traveling wave solutions of the form u(x,t) = <p(x + t-e). However, 
our results amount to the existence of a kind of logarithmic traveling wave behaviour, that 
is a kind of wave solutions that travel linearly if we measure distance in a logarithmic scale. 
This whole issue deserves further investigation. 

1.5. New estimates for the fractional diffusion problem. The study of the sub- and 
super-solutions is strongly determined by the existence of suitable lower parabolic estimates 
for the associated diffusion problem, the Fractional Porous Medium Equation (FPME) 

, m f u,.(x,t) + L s u m (x,t) = for x G M. N and t > 0, 

' \ u(x, 0) = u (x) for x G R N . 

In Section |3j we devote a separate study in the case m > 1 of the behavior of the solution 
when \x\ — > 00, more precisely its rate of decay, for small times t > 0. Our main result says 
that roughly speaking 

u(x,t)~ t\x\^ N+2s ^ 
when |x| is large and t small. The precise result is as follows. 



Theorem 1.4 Let u(x,t) be a solution of Problem (1.10) with initial data uq[x) > such 
that uq(x) > 1 in the ball Bi(0). Then there is a time t\ > and constants C*,R > such 
that 

(1.11) u(x,t) > C*t\x\-( N+2s ^ 

if \x\ > R and < t < t\. 
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The fact that solutions of the FPME with nonnegative initial data become immediately 
positive for all times t > in the whole space has been proved in |12[ 113) . Such result is true 
not only for < s < 1 and m > 1, but also for < s < 1 and m > m c = (N — 2s)+/iV, this 
lower restriction on m aimed at avoiding the possibility of extinction in finite time. 

Precise quantitative estimates of positivity for t > on bounded domains of M. N have been 
obtained in the recent paper [8j. The estimates of that reference are also precise in describing 
the behavior as \x\ — > oo when m < 1 (fast diffusion), but they are not relevant to establish 
the far-field behavior for m > 1. We recall that the limit s — >• 1 with m > 1 fixed we get the 
standard porous medium equation, where positivity at infinity for all nonnegative solutions 
is false due to the property of finite propagation, cf. [21]. This explains that some special 
characteristic of fractional diffusion must play a role if positivity is true. 

We fill the needed gap for some convenient class of initial data that includes continuous 
nontrivial and nonnegative initial data. We give a quantitative version here since it can be 
useful in the applications. 



2 Preliminaries 



2.1 Nonlinear diffusion. The Fractional Porous Medium Equation 

We recall some useful results concerning the porous medium equation with fractional diffusion 
(FPME). We refer to [13] where the authors develop the basic theory for the general problem 

. , f ut = -L s {\u\ m - l u) for x E M. N and t > 0, 

^ ' ' \ u(0,x) = u (x) for x G R N , 

with data uq G L 1 (M. n ) and exponents < s < 1 and m > 0. Existence and uniqueness of 
a weak solution is established for m > m c = (N — 2s)+/N giving rise to an L 1 -contraction 
semigroup. Recently in |14) . it was proved the C 1,a regularity. Positivity of the solution 
for any m > corresponding to non-negative data has been proved in [S]. We give a brief 
discussion on these facts in the Appendix. 



2.2 Barenblatt Solutions of the Fractional Porous Medium Equation 

An important tool that we use in the paper is represented by the so called Barenblatt solutions 
of the FPME. In [23], the second author proves existence, uniqueness and main properties of 
such fundamental solutions of the equation 

(2.2) u t + {-A) s u m = 0, 

taking as initial data a Dirac delta u(x, 0) = MS(x), where M > is the mass of the solution. 
We will give here a short description of these functions and recall their main properties we 
need in the paper. Next, we recall Theorem 1.1 from [23 . 
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Theorem 2.1 For every choice of parameters s £ (0, 1) and m > m c = max{(iV — 2s) /N, 0}, 



and every M > 0, Equation (2.2) admits a unique fundamental solution; it is a nonnegative 
and continuous weak solution for t > and takes the initial data in the sense of Radon 
measures. Such solution has the self-similar form 

(2.3) B M (x,t) = t- a F M (\x\t-P) 

for suitable a and j3 that can be calculated in terms of N and s in a dimensional way, precisely 

N 1 
( 2 - 4 ) — — , (3 



N(m - 1) + 2s ' N(m - 1) + 2s ' 

The profile function FM{r), r > 0, is a bounded and Holder continuous function, it is positive 
everywhere, it is monotone and goes to zero at infinity. 



By Theorem 2.1 there exists a unique self-similar solution B\{x,t) with mass M = 1 of 



Problem (2.2) and moreover, it has the form B\(x,t) = t a Fi(\x\t "). Let BM(x,t) the 



unique self-similar solution of Problem (2.2) with mass M. Such function will be of the form 

(2.5) B M (x, t) = MB\ (x, M m -H) , 
which can be written in terms of the profile F± as 

(2.6) B M (x,t) = M l -( m -^ a r a F ((M 1 m - 1 t)" /? |; 

Moreover, the precise characterization of the profile Fm is given by Theorem 8.1 of |23j . 

Theorem 2.2 For every m > mi = N/(N + 2s) we have the asymptotic estimate 

(2.7) lim F M {r)r N+2s = C x M a , 

r—toc 

where M = f F(x)dx, C% = C\{m, N,s) > and a = (m — mi)(N + 2s) (3. On the other hand, 
for m c < m < mi, there is a constant Coo(m, N, s) such that 

(2.8) limF M (r)r 2s /( 1 - m )=C 00 . 

r— >oo 

The case m = mi has a logarithmic correction. The profile F has the upper bound 

(2.9) F(r) < Cr- N ~ 2s+e , Vr > 
for every e > 0, and the lower bound 

(2.10) F(r) > Cr~ N - 2s logr, for all large r. 

We state now some properties of the profile Fi(r), r > obtained as consequences of formula 



(2.7) that we will use in what follows. Let us consider first the case m > m\. 
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1. Fx attains its maximum when r = i.e. F(r) < F(0), for all r > 0. 

2. There exists K\ > such that 

(2.11) F(r) < K ir -( N+2s \ Vr > 0. 

3. There exists Ki > such that 

(2.12) F{r)>K 2 {l + r N+2s )-\ Vr > 0. 

Similar estimates hold also in the case m c < m < mi, and the corresponding tail behavior 
is different, F(r) ~ r -2s/(i-m)_ This w [\\ have an effect in the different results we get for the 
generalized KPP problem. 

As a consequence, the author also proves that the asymptotic behavior of general solutions 



of Problem (2.1) is represented by such special solutions as described in Theorem 10.1 from 



Theorem 2.3 Let uq = ^ € M+(R N ), let M = [i(R N ) and let Bm be the self-similar Baren- 
blatt solution with mass M . Then we have 

lim \u(x,t) - B M (x,t;M)\ = 

t— >oo 



and the convergence is uniform in 



pV 



2.3 Lower estimates for nonnegative solutions in the case m\ < m < 1 

We use the notations: m c := [N — 2s) /N, $ := l/[2s — N(l — m)] > for m > m c . The results 
we quote are valid for initial data in a weighted space uq G L 1 (M Ar , ipdx), where ip satisfies the 
following conditions: 

Assumption (A). The function £ C 2 {W N ) is a positive real function that is radially 
symmetric and decreasing in |x| > 1. Moreover (p satisfies 

2s 2s 

< ip(x) < \x\~ a for \x\ » 1 and N < a < N H . 

1 — m 1 — m 

We recall now Theorem 4.1 from [8] giving local lower bounds for the solution of the diffusion 
problem. 



Theorem 2.4 (Local lower bounds) Let Ro > 0, m c < m < 1 andletO < uo G L 1 (R N ,ipdx) 
where if is as in Assumption (A). Let u(-,t) G L 1 (K N , ipdx) be a very weak solution to the 



Cauchy Problem (2.1), corresponding to the initial datum uq. Then there exists a time 



(2.13) 



/ fi Di) II,,, II l-m 
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such that 

(2.14) inf u(x,t)>K 1 R Q I ^t^ifO<t<t^ 
and 

\\uo\\ 2s f 

(2.15) inf u(x,t)>Jri °^ Bflo) ift>U. 
The positive constants C m , Ki, K2 depend only on m, s and N > 1. 
The previous estimates, computed for t = t* rewrite as 



(2.16) inf U (x,t)>K 1 Ct m \\uo\\ L i i B Rn )Ro N - 

x ^ B R /2 



Then, if Rq increases, the lower bound will decrease. 

Concerning quantitative lower estimates, we recall Theorem 4.3 from [8]. 

Theorem 2.5 (Global Lower Bounds when m\ < m < 1) Under the conditions of The- 
orem 



2.4 we have in the range m± < m < 1 



C(t) 

(2.17) u(x,t) > ■ |jy +2s when \x\ » 1, 

valid for all < t < T with some bounded function C > that depends on t, T and on the 
data. 

Theorem 2.6 (Global Lower Bounds when m c < m < mi) Under the conditions of The- 
orem \2.4\ we have in the range m c < m < mi 

(2.18) u(x,t ) > C(t)\x\- 2s K 1 - m '> 
if \x\ > R and < t < to. 

The lower estimates for exponents m > 1 need a new analysis that we supply in the next 
section. 



3 Lower parabolic estimate in the case m > 1 

We consider the FPME equation (with no reaction term) for x G H N and t > with nonneg 
ative and integrable initial data 

(3.1) u(x,0) = uq(x) , 
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and we also assume that uq is bounded and has compact support or decays rapidly as \x\ — >• oo. 
We want to describe the behavior of the solution u(x,t) > as |x| — > oo, more precisely its 
rate of decay, for small times t > 0. We take m > 1 since the study of positivity for rn < 1 
was dealt with in previous results. 



The first step in our asymptotic positivity analysis of solutions of (2.2)-(3.1) is to ensure that 
solutions with positive data remain positive in some region. We only need a special case that 
we quote next, based on the positivity results of |8j. 



Theorem 3.1 (Local lower bound) Let u be a weak solution to Equation (2.2), corre- 
sponding to uo £ L 1 (M d ). Then there exists a time 

(3.2) U := C R 2s+N ( m ~V \\u4l^ , := l/[2s + N(m - 1)] > 0, 

such that for every t > t* we have the lower bound 



(3.3) 



xeB B/9 - t N # 



inf u{x, t) > K 



valid for all R > 0. The positive constants C and K depend only on m, s and N , and not on 
R. 



We may now state and prove the main lower estimate with precise tail behaviour, which is 
based on a delicate subsolution construction. 



Theorem 3.2 Let u(x,t) be a solution with initial data uo(x) > such that uq(x) > 1 in the 
ball Bi(0). Then there is a time t± > and constants C*, R > such that 

(3.4) u(x,t) > C*t\x\- (N+2s 1 

if \x\ > R and < t < t\. 



Proof. We consider the FPME for m > 1 and initial data uq is 1 in the ball of radius 2. We 
use the previous result to prove that in the ball of radius 1/2, then u{x,t) > 1/2 for < t < to, 
a time that is calculated from the formulas above. 

• We want to construct a sub-solution of the form 

U m (x,t) = G(\x\) +t m F m (\x\). 

We want to choose G > and F > in such a way that U will be a formal sub-solution of the 
FPME in a domain of the form Q = {\x\ > 1/2, < t < h}, i.e., we want U t + (-A) s U m < 
in Q. Note that 

U t = (G(x) + t m F m {x)) {1/m) - 1 t m - 1 F m (x) < F(x). 
We also have, with L s = (— A) s , 

L s U m = L s G{x) + t m L s F m {x). 
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We take F positive, smooth and F{r) ~ r ~(N+2s) as r qq to get the desired conclusion after 
the comparison argument: u(x,t) > U(x,t) > ctr~^ N+2s ' > if r is large and t ~ 0. For later use, 
let us say that F < C 2 r~ {N+2s) for r > 1/2. Since m > 1 we can choose F smooth so that 
L s F m = 0(r~( N+2s ' ) ) for r > 1/2 (use the asymptotic estimates like the first lemma in [HJ) 

We will take G(r) = for r = \x\ > 1/2 so that U(x,t) = tF(x) there. If G is also smooth 
we have L S G bounded and L S G ~ —C\r~^ N+2 ^ as r — > oo. By contracting G in space, 
G(x) = G(kx), k > 0, we may then say that L S G < — Cir~( N+2 ^ for r > 1/2. Then we will 
have for r > 1/2 

U t + L s U m <F + L S G + t m L s F m < C 2 r- {N+2s) - C 1 r"( iV+2s ) + t m L s F m < 

(C 2 + £ )r-^+ 2s ) - C 1 r"( 7V+2s ) < 
if C\ > C*2. We can choose G large so that C\ is large enough. 

• We now want to use the viscosity method to compare u(x, t) with U(x, t) in the region 
Q = > 1/2,0 < t < ti}, and this will prove that U(x,t) < u(x,t) in Q. Apart from the 
sub-solution condition that we have checked, we need suitable comparison of the boundary 
conditions at r = 1/2 

U(l/2,t) = tF{\/2) < 1/2. 

This ends the construction if the comparison result is justified. The contradiction argument 
at the first point of contact between u and U will be justified as in [8] (where it was applied 
to fast diffusion equations of fractional diffusion type) if the solution we have is a bit smooth: 
ut and L s u m must be continuous and the equation must be satisfied pointwise there. This 
regularity is true and the proofs are under study now. 

Alternatively, we may use Implicit Time Discretization with a sequence of approximations. 
The justification of the method in the elliptic case is done in the paper in collaboration with 
Volzone |25j on symmetrization techniques. □ 

Remark. The level uq(x) > 1 in the ball -Bi(O) can be replaced by uq(x) > e > in any 
other ball by means of translation and scaling. In this way the result is true for all continuous 
and nonnegative initial data uq, of course nontrivial. 



4 Evolution of level sets of solutions to Problem (1.1) 



In this section we start the proof of the main result of the paper on evolution of level sets with 
exponential speed of propagation. In a first step we prove the convergence to zero on outer 
sets. Since the decay assumption on the initial data is the same for mi < m < 1 and m > 1, 
we will threat both cases, as well as m = 1, in the following lemma. 



Lemma 4.1 We consider m > m\ and let u be the solution of Problem (1.1) with initial 
datum uq(x) € L 1 (R N ) 



(4.1) 



< uq < 1. We assume that uq satisfies the decay property 
uo(x) < C\x\-( N+2s) for all x £ R N . 
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Then, for a > 03 if m > 1 (respectively, for a > 02 if m\ < m <1), we have 
(4.2) u(x,i)->0 as t -> 00 

uniformly for \x\> e at . 



Proof. We consider the solution u(x, t) of the linearized problem 

u t + L s u m = f'{0)u, u(0,x) =u Q (x). 
Since / is a concave function, we have f'(0)s > f(s) Vs G [0, 1], and thus u is a super-solution 



of Problem (1.1), which implies the upper estimate 

u(x,t) < u(x,t), for t > 0, x G R N . 



Next, we define v(x, r) by 
(4.3) 

and new time 

(4.4) r = 

(4.5) r = 



t>(x,r) = e _/ ' (0)i u(x,i), 



(m-l)/'(0) 



1 



3 (m-l)/'(0)t 



1 



-(l-m)/'(0)t 



if m > 1, 



if m < 1, 



(l-m)/'(0) L 

and t = t for m = 1. It is immediate to check that r) is a solution of the FPME 



(1.10) with initial datum vq = Uq. Let Bm(x,t) the Barenblatt solution with mass M of the 



FPME, as defined in Section 2.2 By virtue of the properties of the Barenblatt solutions and 



assumption (1.5) on the initial data, we conclude there exists M > big enough and tq > 
such that 

v (x) < B m (x,t ). 



By the Maximum Principle 

v(x,t) < B m {x,t + t ), Vx G 



t > 0. 



Now, using the characterization of the decay of the Barenblatt profile given by relation (2.7), 



we obtain that there exists K\ > such that F(r) < K\r ( Ar + 2s ) ) for all r > 0. We obtain the 
:r estir 

u(x, t) 



following upper estimate on the solution u of Problem (1.1): 

= ef'^ t v(x,T) 
< e f ' i0)t B M {x,T + T 



e /'(0)* (r + ro) - Qj p ( | x | (r + ro) -/3 ) 



< ef'^ir + t )-« Ki (\x\(t + ro)^))-^) 
= K 1 ef'^ t (r + T ) 2 ^\x\-( N+2s \ 
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Case m > 1. In order to continue the estimate, we remark that for large times t, the term 
T 2/3s an influence on the result only in the case m > 1. Then (r + tq) 2 ^ s < e^" 1-1 - 1 ^ for 
large £. Let us assume that |x| > e' 7 '. Then one has 

u(x,t) < CK 1 e f ' ( - 0)t T 2l3s e-^ N+2s)t = CKie [/ ' (0,+2/ ' (0)(m " 1),3s " ff( ' v+2s),i . 

We want to have /'(0) + 2/'(0)(m - l)/3s - c^A 7 " + 2s) < 0, which is just the condition 

l + 2(m-l)/3s 
> N + 2s f (0) = 

We have obtained the convergence of u(x,t) to as t — > oo, for > e CTi . 

Case m < 1. In this case, the term (r + tq) 2i3s is bounded for every t > as we can see from 



(4.5). As before, we assume \x > e ut . Then, we get 

u(x,t) < CK l ef'^ t e-^ N+2 ^ t = CK 1 e^'^-^ N+2 ^ t . 

For a > 02 = > the exponent is negative /'(0) — o"(A^ + 2s) < and we obtain the 
convergence of u(x, t) to as t — > oo. 

□ 



Lemma 4.2 VFe consider m c < m < mi. Let u be the solution of problem (1.1) wzi/i initial 
datum uo(x) € L 1 (IR Ar ), < uq < 1 and we assume uq satisfies the decay property 

uq{x) < C\x\~ 2s ^~ m ^ for all x G R N . 

Then, for a > o\ we have 

u(x, t) — > 0, i — >• oo 

uniformly for \x\ > e 17 *. 



Proof. The proof follows the same as in Lemma 4.1 since the Barenblatt solution Bj\j(x,t 
T~ a F(\x\r~P) of the diffusion problem satisfies F(r) ~ j-~ 2s /( 1_m ) according to Theorem 
Therefore, we obtain the estimate 

u(x, t) < e f '^\T + r )- a Ki(|x|(r + ro )-^))-2V(l-m) 

= ^ ie /'(0)* (r + ro) -l/(l- m )| xr 2 S /(l- m)> 

Since m < 1, the term (r + To) -1 ^ 1- m ) is bounded and then, for \x\ > e°"* we obtain 

u(x,t) <^ ie /'(0)t-2.«rt/(l-m)_ 

For o~ > T^r/'(0) = o"i we obtain the desired convergence to as i — > oo. 



2.2 



□ 
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Remarks 



I. When m = 1 we recover the minimal speed a 2 = f'(0)/(N + 2s) obtained by Cabre and 
Roquejoffre in [TT]. The proof is similar, but in the nonlinear case we have to make an 
exponential change of time variable. Note also that we only use the decay properties of the 
fundamental solution. 

II. The value of the critical exponent 02 can be easily obtained from the following formal 
study of the level lines of u(x, t). Thus, the set {u(x, t) ~ e} can be written in terms of v(x, r) 



defined in (4.3) as 
(4.6) 



-/(°>*t;(x,r) 



By Theorem 2.3, v(x,t) behaves like the Barenblatt solution of the Fractional Porous Medium 
Equation (2.1 ): 



v(x, t) ~ B(x, t) = t a F(r), F(r) ~ r 



-(N+2s) 



XT 



From [23], we know that B(x,t) ~ T -^(N+2s)^-(N+2s) ^ thug _ T 2p s \ x \-(N+2s) _ At 

this moment, Co) implies e /,(0) V 2/3s |x|" (Ar+2s) ~ e. 



For instance in the m > 1 case, it follows that 



_ ^ e (l+2/3s(m-l))/'(0)t 



l/(JV+2s) 



l + 2/3a(m — 1) cl / n \, 
e N + 2s ( ' 



1 _|- 2/3s(m — 1) 

and we deduce an exponential behavior of the level sets \x\ ~ e' 73 *, where 03 = Ar i o ^ /'(0). 



Similarly, in the mi < m < 1 case, we get that 



N + 2s 



.r I -e 



/'(0)t 



l/(JV+2s) 



f(0) 
iV + 2s' 



5 Evolution of level sets II. Convergence to 1 on inner sets 



In this section, we will prove the convergence to 1 of the solution u(x, t) of Problem ( 1.1 ), i. e., 
the second part of the statements of our main theorems |1.1[ |1.3[ and |1.3| 



5.1 Case m > mi 

We will present this case in full detail. The proof for the case m c < m < mi being similar, we 



will sketch it at the end of this section. We have iV > 1, s E (0, 1), m > mi, f satisfies (1.2), 
and o"2 



f(Q) 

N+2s 



as defined in (1.6). 



Proposition 5.1 Let N > 1, s G (0,1), mi < m, f satisfying (1.2). Let u be a solution of 



Problem ( 1.1 ) with initial datum < Uo(-) <1 } Uq^= 0. Then for every a £ (0, 02), u(x, t) — > 1 
uniformly on {\x\ < e at } as t — > 00. 
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Proof. We fix a £ (0,02). Proving the converge of u(x,t) to 1 is equivalent to proving the 
convergence of 1 — u m to 0. Therefore, we fix A > and we need to find a time t\ large enough 
such that 1 — u m (x, t) < A for all t >t\ and \x\ < e at . 

• Let us accept for the moment the following lower estimate that will be proved later as 
Lemma 5.4 given v G (cr, 172)1 there exist e S (0, 1) and to > such that 

(5.1) u > e, for t > to and \x\ < e vt . 

We now proceed with the last part of the argument, where the effect of the nonlinear diffusion 
is most clearly noticed. We take t\ > to and consider the inner sets where 

e < u < 1 for (x, t) E SI/ := {t > h, \x\ < e ut }. 

Let v = 1 — u m . Then v satisfies the equation 



1 



(5.2) -(l-v)^- 1 v t + L s v + f{u) = 0, 

m 

that we write in the form 

(5.3) a(x,t)v t + L s v + b(x,t)v = 0, a(x, t) = — -u 1_m , b[x,t) = 

m v 

Moreover, we estimate a(x, t) as follows 

ao = — e 1_m < a(x,t) < a\ := — in SI/, if m < 1, 
m m 

respectively, 

ao = — < a(x, i) < a\ := — g 1- " 1 i n ^7, if m > 1. 
m to 

We argue similarly for b(x,t) in $7/: 



1 — u m (1 — u)mt, m 1 



where 

60 = if m < 1 and 6 = — if m > 1. 

ml — e ml — e 

In particular, v satisfies 

(5.4) a(x, t)vt + L s v + 6qv < in f2/. 



• We look for a super-solution w to Problem (5.3) that will be found as a solution to a linear 
problem with constant coefficients, and we also need that wt < 0. More precisely, we consider 
w solution of the concrete problem 



(5.5) 



a\Wt{x, t) + L s w(x, t) + bow 
w(x, t\) = 1 + c^M 7 







for x £ 
for x £ 



and t > 0, 
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where the exponent 7 taken such that 
(5.6) 



0< 7 :=- 6 ° <2s. 
v a\ 



We can eventual consider a smaller e for this inequality to hold. Equation (5.5) is linear, the 
solution can be computed explicitly 



w{x,t) = e °i w{x,t), 
where w(x, r) solves the linear problem 

w T (x, t) + L s w(x, t) = 0, w(0) = 1 
We observe that w can be written in the following form 



1 



t = — (t-ti), 
ai 



(5.7) 
where 



w(x, t) = 1 + -^-U(x, T + n), 

^2 



01 = 2? A = 2? 



U(x,t) = T ai F(\x\r~ Pl ), 

is the self-similar solution of the linear problem 

U t (x,t) + L s U(x,t) = 0, U(x,0) = |x| 7 . 

The properties of the self-similar solutions t/(x, r) deserve a separate study, which is done in 
detail in Section |6j Thus, by Lemma 6.2 the profile F is non-decreasing and U (x, r) has a 
spatial decay as |x| 7 for large |x|t -1 / 2s : 



(5.8) 



C 2 |x| 7 < U(x, t) < Ci|x| 7 for all |x|T _1 / 2a > K x . 



We will consider a suitable delay time t\ in the definition of w stated in (5.7). In what follows 
we will use the notation 77 = \x\t~ ? x . We check that the derivative wt is negative: 



w t (x,t) 



d_ 
dt 



e % ( - t - tl \l + C% 1 U{x,T + Ti)) 



e Q i 



•^(1 + C^\r + nrF(rj)) + C^(t + n)" 1 " 1 (aiFfa) - Pi V F'(r))) ^ 

= e"^^"* ^- [-6 C 2 + (-6o(r + n) + ai)(r + n) 01-1 ^) - ^(r + n)" 1 " V'M] 
aiU2 

Since ^'(77) > for all 77 > 0, we get that wt(x, t) < for all t > t\ if r + n > ai/bo, which is 
true for a suitable choice of t\. 



Now we can compare w and v by applying the Maximum Principle stated in Lemma 8.1 



of the Appendix , as in [TT] . Define v = v — w and ensure the hypothesis of the Lemma are 
satisfied. 
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(HI) We check that w(x,h) > v(x,h) for all x G R N : 

w(x,ti) > 1 > v = 1 - u m , Vj; G 



(H2) We check that u> > w in ((ti, oo) 



\ SI/, that is f > t\ and |x| > e . At this point, 



we use the estimates (5.8). We ensure that e ut > K\(t + t\) 1 / 2s for all £ > t\, which is true 
by choosing eventually a larger t\. Therefore 



w 



{x,t) = e^ {t ~ tl] w{T,x) > e~£ ( ' _tl) (l + -^C 2 \x\ r ) 

(^2 



> e 



1 + e 7i ") > 1 > v(x,t) for all t > h, \x\ > e 



since 7 satisfies (5.6). By the previous computation v < in ((ti,oo) 



(H3) Next step is to prove that v is a sub-solution of Problem (5.4). Indeed, we have that 
a(x, t)vt + L s v + bov = a(x, t)v t + L s v + b$v — [a\Wt + L s w + bow] + {a\ — a(x, t))wt < in 
By Lemma 8.1 we obtain that v < in [ii,oo) x M. N for t% taken to be large enough. Thus, 

Ci 



v(x, t) < w(x, t) 



e ^ {t ' tl] (I + C^U(x,t)) < e' b ^ {t ~ tl] (I + ^-\xP ■ 



+ ^ e ^c7 < A 



Let us consider the inner set (x,t) £{t> t x , \x\ < C x e ut }. We have 

ft' 

for Ca small enough and t\ large enough. 

Finally, since a < v then e at < C\e vt for every t > t\ with large enough, and the previous 
inequality implies that 

1 - u m {x, t) = v(x, t) < A for t > t x , \x\ < e at , 

which concluded the proof of the uniform convergence to the level u = 1. □ 

To complete the proof of the result of this subsection, we need to supply the proof of the 
lower estimate (5.1). This will be done in three steps, 

Step I. Starting with arbitrary initial datum < no < 1, uq ^ 0, we obtain a lower bou nd f or 
u with the desired tail u > c \x\~( N+2s ^ for large \x\. The result corresponds to Lemma 



5.1 



Step II. We prove that given an initial data taking the value e in the ball of radius po an d 
decaying like that |a;| _ ^ +2s ) for large the corresponding solution of Problem (1.1) will be 
raised to at least the same level e in a larger ball p\ and in a later time that is estimated. The 



sizes are important. This will be Lemma 5.2 



Step III. By combining the previous two results, we conclude that u > e on the inner sets, for 
a certain e > 0. This will be Lemma 15.31 and Lemma 15.41 



Steps II and III follow the ideas of [llj in the linear case, with a long technical adaptation to 
nonlinear diffusion. 
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Lemma 5.1 (Long Tail Behaviour) Let N > 1, s G (0,1), m > mi, f satisfying (1.2) 



and a G (0,02). Let u be the solution of Problem (1.1) with initial datum u(0, •) = uq, where 



< Uq < 1, Uq 7^ 0. T/ien /or any /ixec? to > i/iere exist e G (0, 1), ao > 0, po > 1 such that 



u(x, t) > vq(x) :- 



a \x\-( N+2s \ 

-(N+2s) 

e = a p y 



\x\ > po, 
\x\ < po, 



for all t G [to, 2to], x G 



N 



Proof. We recall that o~i = f'(0)/(N + 2s). The idea is to view u the solution of Problem 



(1.1) as a super-solution of the homogeneous problem with the same initial datum uo, that is 



the FPME. Therefore, 

u(x,t Q + t) > u{x,t), Vt > 0, x G R N , 
where u is the solution of the FPME with initial datum uq 



(5.9) 



u t (x, t) + L s u m (x, t) = for x G 1^ and t > 0, 
u(x, 0) = uq(x) for x G M^. 



We will estimate u from below by using the local and global estimates on the FPME given in 



Theorem 2.4 and Theorem 2.5 for m < 1, respectively Theorem 3.1 for m > 1. The decay in 
case m = 1 is well known, see Section [6] for a review. More exactly, in all cases m > mi, there 
exist a time T > and constant R > such that 

u(x,t) > C(t)|zr (JV+2s) , V|x| > R, 0<t<T. 

Then, for a fixed G (0,T) which also satisfies i* < io ; we can find a Barenblatt solution 
BM(x,t) and a time t 2 > such that 

u(x, i*) > S M (x, t 2 ), Vx G R N , 

and therefore, by the Comparison Principle 

u(x,t + Q > B M (x,t + t 2 ), Vx G M^, t > 0. 

In particular, we can choose e > such that 



u(x, t) > Vo(x) :- 
for all x G R N , t G [i ,2i ]- 



ao|x 



-(iV+2s) 



-(N+2s) I I ^ 

aoP , \x\ < Po, 



□ 
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Lemma 5.2 (Positivity for a sequence of times) Let m > m\. For every a < o~2 there 
exist to > 1 and < eo < 1 depending only on N, s, f and a for which the following holds: 
given po > I and < e < eo, let ao > be defined by aop ^ +2s ^ = e, if we take 



(5.10) v (x) 



ao\x\~( N+2s \ \x\ > po, 

-(N+2s) I I . 

e = aoPo > \ x \ ^ Po, 



then the solution v of Problem ( |1.1[ ) with initial condition vq satisfies 
(5.11) v(x,kt ) > e for \x\ < p e akt °, 

for all k G {0,1,2,3,...}. 



Proof. Case m > 1. 

I. Preliminary choices. From the beginning we fix a G (0, 02). We will do a very detailed 
analysis of the case k = 1, which is then iterated for the rest of values of k. We take 5 G (0, 1) 
small enough such that 

For example, take 5 such that 

m 1 



(N + 2s)S 2 



o"2 + max{<7, 7V(m — l)/3o"2}) • 



This choice will be explained later. Next we take to sufficiently large depending only on N, s, uo 
and a such that 

(5.13) e^* (1 + (t /t 1 (1)))- n ^ > K 3 , (K 2 /2K 1 ) 1 ^ N+2 ^ e ^™ to > e at ° . 

where n(l) = e 1_m C2 with C2 a positive constant that we state explicitly later, and K 2 < 2K\ 
are constants describing the properties of the profile F\ of the Barenblatt function with mass 
1 given in (2.11) and (2.12), and we recall for convenience that 

K 2 (l + r^ 23 )- 1 < F{r) < K ir -^ N+2s \ Vr > 0. 

Define now eo G (0, a) by 

eo = <5e-/'(°)'o. 
Clearly, eo < S. Now, we fix < e < eo and po > 1. 



II. Construction of sub-solutions to Problem (1.1). Let w be a solution of the problem with 
linearized reaction 



(5.14) 



w t {x, t) + L s w m (x, t) = ^f-w for x G 1^ and t > 0, 
w(0, x) = vq(x) for x G H N . 
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We define w(x,t) by 

with a new time 
(5.15) 



w(x, t) = e s w(x, t) 



(m-l)/'(0) 



,(m-l)f'(0)t 



1 



if m > 1, 



so that r = t in the limit m = 1. Then, w is a solution of the Fractional Porous Medium 
Equation with initial datum vq 



(5.16) 



w T (x, t) + L s w m (x, t) = for x E 1^ and r > 0, 
W(x, 0) = uo(^) for x E M^. 



III. Comparison with a Barenblatt solution. Lower bound for v(x,to). We prove that there 
exist Mi > and t\ > such that 



(5.17) 



vq(x) > B Ml (x,n), Vx E 



where Bmx(x,t) is the Barenblatt solution of Problem (1.10) with mass Mi given by Theorem 

OLD 

(5.18) B Ml (x, t) = M 1 B 1 (x, Mf-V) . 
Now, Bm 1 {x,t) can be written in terms of the profile F\ as 

(5.19) B Mi (x,t) = M 1 1 - {m ~ l)a T- a F 1 ((Mf-V)"^ |. 

We will use the properties of the profile F\ stated in ( |2.11[ ) and ( |2.12 ). With this information, 
we will find the constants M\ > and t\ > such that inequality (5.17) at the initial 
time holds true. For \x\ < po we have that Bm 1 (x,ti) < M\ ^ m 1 ^ a T 1 _a Fi(0). Note that 
1 — (m — l)a = 2 (3s > 0. We impose the first condition 

(5.20) M 1 2/3s r 1 - a F 1 (0) < e. 
Let \x\ > po- Then, 

B M Ax,r) < Mf s r^K x ([M?-\yP |^Q- (iV+2s) 
In order to use this the inequality for large \x\ we also impose the condition 

(5.21) K x M{ +2p * m - x) T$ f>a < oq, where a = ep£ +2s . 



Conditions (5.20) and (5.21 ) are sufficient for inequality (5.17) to hold. Under such restrictions 
the st 

(5.22) 



the stated inequality (5.17) holds true. Then, by the Comparison Principle we get 

Bmi (x, t + ti) < w(x, t), for all \x\ E M. N , 



t > 0. 
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Putting equality in the inequalities (5.20) and (5.21) we get 
(5.23) Mi = Cl ep£, n = c 2 e l - m pl s , 

(with c\ , C2 positive constants not depending on e or p). We can easily see that the expres- 
sions are dimensionally correct. The constants are given by c\ = Fi(0)~ 2s ^ N+2s ^ K 1 N /( N+2s ) j 

C2 = Fl ^(i+2(m-i)M/(W+2s)) K - 2s/ ( N+2s \ In particular, (M^nf = c m , with c 3 = 
Jm-l)flJ 



1 



'2- 



Since v o < e in R then u;(x, r) < e for all x G M , r > 0, and then in terms of w(x, t) we 
obtain the following bound 



< w(x, t) < ex'Of < e /,(0)<0 e = 5, Vt < t . 



Since ^^-w < f(w) for 0<w<5,wis& sub-solution of Problem ( 1.1 ) in [0, to] x ^ N ■ By the 
Comparison Principle and estimate (5.22) we obtain that at the moment to 

(5.24) 



/(<5)j /(<5)j- AT 

v(;t ) > w(-,t ) = e^ to w(-,T ) > e^r t0 B Ml (-,T + n) in TOA 



where we use the notation tq = r(to) defined by (5.15). 

IV. We will now prove that estimate ( |5.24 ) with the choices (5.23) for M\ and T\ implies the 
lower bound stated in Lemma 5.2 in the case k = 1, m > 1. Indeed, we have (*) 



f(6) , 

v(x,t ) > e s t0 B Ml (x,T + n) 

= e^M 2 / s (r + n)^ fM-^-^l 



to + n, 



/(£) 



> e^M^'fo + n)- a K 2 (l + (M- (m - 1)/3 (r + 7 1 r /3 l^< x +~ >M 

^e^Mf'fro + n)-" f 1 + C3 (7V+2S) (1 + (ro/n))^^ (M/p ) (7V+ ^ 



r -(m-l)/3. 



Our aim now is to be able to continue this estimate until we reach a bound v\{x) of the form 
(5.10) for the same e and a different parameter p\. We will choose some p\ and then check 
that the lower bound for v(x, to) is larger than e at |x| = pi- In order to simplify the estimate 
of the last parenthesis in formula (*), we will impose the condition 

cs {N+2S) (1 + (Wn))"^ (Pi/P0) {N+2S) > 1 , 
and then we only need to have 

(5.25) v( Pl , to) > (K 2 /2)e^Ml +2 ^ s (ro + n) 2 ^"^ > * 

The first condition is equivalent to 

(Pi/ P o) iN+2s) > cf +2s (i + (Vn)^ 
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while, taking into account that M l 
that 



l+2(m-l)/3s 2/3s 



r^ ps = ao/K\ and ao = ep^ +2s , the second means 



(5.26) (pi/po) (N+2s) < {K 2 /2K{)e^ (1 + (r /n)) 2s/3 . 

Both conditions are compatible iff 



e^*o (1 + (ro/n))"^ > 



2Ki 



,N+2s 



Now recall that n depends on po by (5.23), and t\ is bounded below by ti(1), the value for 



p = 1 . We see this condition as a way of choosing to . Using the fact that for large t we have 

r~e( m - 1 ^*/(m-l)/'(0), 
we easily see that for large to the left-hand side looks like 

Ce ^t -N(m-l)l3f>(0)t ^ 

hence, the compatibility condition can be solved if f(5)/5 > N(m — l)/3/'(0). Since 5 is small 
enough so that f(5)/5 ~ f'(0), this means that we need N(m — l)f3 < 1 which is true. We 
conclude that there exists to > large enough such that 



(5.27) 



e^'°(l + (To/Ti(l))) 



-N/3 



This choice of to is independent of pq. 



Once this is guaranteed, we choose the largest possible pi satisfying (5.26), which is 
(5.28) 



^ = (K 2 /2K 1 ) 1 ^ N+2 ^ e TO^° (i + {TQ / Tl) fsp/(N + 2s) . = L 
Po 



V. With this choice of pi and to, estimate (5.25) holds. Going back to Point IV above, we 
have 



v(x,t ) > e"« to B Ml (x,T(t ) +n) > e, for |x| = pi 
and thus, since the profile .Fi is non-increasing we get that 

v(x, to) > e, V|x| < p\ . 



Using (5.28) we get that 



v{x,t )>ep^ +2s \x\~ {N+2s \ y\x\>Pi. 
Finally, we define a\ := ep^ +2s and thus v(-, to) > vi(-) where vi is given by the expression 



vi(x) 



ai|x|-^ +2s ), \x\>pi; 

~(N+2s) | | , 

e = a 1 p 1 \x\ < pi. 
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The proof is complete for m > 1 and k = 1 (see Figures [2] and [3] for the construction of v±). 

VI. The iteration. We are now ready to address the next delicate step. Once we have proved 
that v(x, to) > vi(x) for all x E M , where «i is defined above, we apply the same proof and 
result to obtain 

v(x, 2to) > (solution of KPP with initial data vi(x))(to) > V2(x), 

where v^ipc) has the same construction as vq(x) and v\(x) but with parameters pi and 0,2- Since 
pi > po the previous choice of to is still valid to get to a similar conclusion. The argument 
continues for all k = 3,4, .... 

Let us check more closely the quantitative part of the iteration in order to get an improvement. 
In the process we keep e fixed but me replace po by p^, k > 1, so that the formula (5.28) 
becomes 

^±1 = Lk := (i^A^) 1 /^ e(^ t0 (1 + (Wny^ • 
Pk 

As k — ^ 00 we have p& —¥ 00, hence T\{pk) — > 00, and the last quantity tends to 

Loo = {K2/2K 1 ) l ^ N+2s ^ e W^ t0 . 
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Finally, if we are given some a < o<i = f'(0)/(N + 2s) we can change the definition of to so 
that we also have > e at °. The conditions we put on 5 and to can be summarized in (5.12) 
and (5.27), and they are independent on the parameter tj-, Pk of the iteration. This ends the 
proof for m > 1. 

Case m < 1. The outline of the proof is similar to the case m > 1. We explain the differences 
that appear. The new time r is introduced via 



(5.29) 



1 

(I - m)f'(0) 



1 



-(l-m)/'(0)t 



if m > 1. 

= !/((! — m) /'(())). This property 



Therefore, for each t we have a new bounded time r(t) < r D 
allows us to simplify the choice of to as follows: condition (5.27) is satisfied if 

(5.30) 

where t 00 /t(1) = e 

Summing up: we take 5 small enough such that 

/(*) 



>if 3 (l + (Wri(l))r, 
" (02(1 — m)/'(0)) _1 is a constant independent of po,pi- 



(N + 2s)5 



> a 



and to such that 

(5.31) e^ t0 >K 3 (1 + (WTia)))^ . 

The rest is essentially the same. 



(K 2 /2Ki) 1/( - N+2s) eW^&js to > e <rto_ 



□ 



Lemma 5.3 ( Expansion of uniform positivity for all times) Let N > 1, s G (0, 1) ; 

Let to > /rom Lemma |5.2| 



mi < m, f satisfying (1.2) and cr G (0,0-2). Let to > /rom Lemma |5.2| T/ien /or every 
measurable initial datum uq with < uq < 1, uo ^ 0, there exist e G (0, 1) and 6 > ^ot/i 
depending on uq) such that the solution u of Problem (1.1) with initial datum u(0, ■) = uq 
satisfies 



u(x, t) > e for all t > to and \x\ < be 



at 



Proof. Let to defined in Lemma 5.2 Then by Lemma 5.1 there exist e>0, ao>0, po>l 
such that u(x, t) is bounded from below by a function vq with the long tail behavior at infinity 



u(x, t) > vq(x) :- 



ao\x\ 



-(N+2s) 



\x\ > po, 



-(N+2s) I 1 , 

aoPo , \x\ < po, 



5.2 



for all x G R N , t G [t , 2t ]. In this way vq can be taken as the initial datum (|5.10l) in Lemma 
We make ao smaller, if necessary, to have that e = aop ( N+2s ^ < eQ; where eo is given in 



Lemma 15.21 
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Therefore, by applying Lemma |5.2[ the solution u will be raised an e at a large time tq + to 



and this holds true for all To G [to,2io]. More exactly, by (5.11), for every k = 0, 1,2, ... one 
has 

u(x,t + kt ) > e for all \x\ < e akt °p , r G [to,2t ] 

which rewrites as 

(5.32) u(x,t) > e for all \x\ < e akt °p , t G [(k + l)t , (fc + 2)t ]. 

But for i G p + l)i , (fc + 2)t ] we get e CTfc *° = e CTfc *o-^ e ot > e -2<7t 0e ^ and then ( | 5 32 [ ) implies, 
in particular, that 

u(x,t) > e, for all \x\ < e' 2cjto e at po, t G [{k + l)t , (k + 2)t ]. 

Since the union the intervals [{k + \)tQ,(k + 2)to] with k = 0,1,2,.... cover all [to,oo), we 
deduce that 

u(x, t)>e if t > t and \x\ < poe~ a2to e r7t . 
The proof of the lemma follows by denoting b = poe~ a2t ° . □ 



. Let u be a solution 

of Problem (1.1) with initial datum < uq(-) < 1, uq 7^ 0. TTien for every a < 02 there exist 



Lemma 5.4 Let N > 1, s G (0,1), / satisfying (1.2). Lei 02 



go) 



e G (0, 1) and t>0 such that 

u(x, t) > e /or all t >t and \x\ < e CT *. 



Proof. We apply Lemma 5.3 with cr replace by o 7 £ (<j, (T2). Since e at < be"' 1 for t large, where 



b is the constant in the statement of Lemma 5.3 we deduce that 



u(x, t) > e for t > t and \x\ < e a . 



□ 



5.2 Case m c < m < mi 

In a similar way, we can prove the convergence to 1 on the inner sets also in the range of 
parameters m c < m < mi. 



Proposition 5.2 Let N > 1, s G (0,1), m c < m < mi, / satisfying (1.2). Let o~i = 
l ~ m -f'(0). Let u be a solution of Problem (1.1) with initial datum < uq(-) < 1, ^ 0. 



2s 



T/ien /or every a G (0, o"i), n(x, t) — )• 1 uniformly on \\x\ < e } as t — > 00. 



Proof. We argue in a similar way as in the case m > mi proved in Proposition 5.1 The 
difference appears when obtaining the positivity on inner sets. To this aim, we start with 



nontrivial initial data < uq < 1 and we prove the analogue of Lemma 5.3 The key ingredient 
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is to use the quantitative lower estimates for the solution u(x, t) Fractional Fast Diffusion 
Equation stated in Theorem |2.6| to obtain an estimate of the form 



u(x,t) > v (x), Vi G [t ,2t ], x G R N , 

where vq(x) is defined as 



(5.33) v (x) 



aokr 2s/(1 ~ m) , \x\>po, 

-2s/(l-m) I I ^ 

e = aoP i Fl < Po , 



Afterwards, we can prove an analogue result to Lemma 4.1 starting with initial data of the 



form (5.33). Since the Barenblatt solution has a long tail decay of the form \x\ 2s K l m \ then 



we find M\ > and t\ > such that 

vo(x) > B Ml (x,n ), VxGM^. 

□ 

6 The linear diffusion problem 

We will need a number of facts about the linear diffusion equation for < s < 1, 
(6.1) U t + (-A) S U = iovx£R N and t > 0. 

This problem has been studied, mainly by probabilists ([2J[6]), see also [22], and many results 
are known. When considering initial data U$ G L 1 (M Ar ), or more general, 



(6.2) U(0,x) = U o (x) forxG 



the solution of Problem (6.1)-(6.2) has the integral representation 



U(x,t)= / K s (x - z,t)U (z)dz , 
Jr n 

where the kernel K s has Fourier transform K s (^,t) = e - ^' 2 "*. If s = 1, the function K\{x,t) 
is the Gaussian heat kernel. 

6.1 The fundamental solution. Further results on the asymptotics for large 

|.r| 

We need some detailed information on the behaviour of the kernel K s (x, t) for < s < 1. In 
the particular case s = 1/2, the kernel is explicit, given by the formula 

K 1/2 (x,t) = C N t(\x\ 2 + t 2 )- {N+1)/2 - 
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In general, we know that the kernel K s (x,t) is the fundamental solution of Problem (6.1), 
that is K s (x,t) solves the problem with initial data the Delta function 

lira K s (x,t) = 5(x). 
t— >o 

It is known that the kernel K s has the form 

K s { X ,t) = t- N ' 2s f{\x\t- l l 2s ) 

for some profile function, /(r), that is positive and decreasing, and behaves at infinity like 
f{r) ~ r -( N+2s \ cf. 0. 

We perform now a further analysis of the properties of the fundamental solution. Our aim is 
to prove the following result. 



Proposition 6.1 For every s E (0,1), the fundamental solution K s {x,t) of Problem (6.1) is 
a increasing function in time 

— K s (x,t) > 0, for all large values of \x\/t l ^ 2s . 

This property is known to be satisfied for the fundamental solution of various types of diffu- 
sion equations of evolution type: the Gaussian profile for the Heat Equation, the Barenblatt 
solution for the Fast Diffusion Equation. 

The analysis of the derivative 4rK s (x,t) involves not only the characterization of the profile 
/ for large r, but also a similar property for the derivative /'. In fact, we will prove that 
f(r) and rf'(r) have the same behavior for large arguments. This is due to the power decay 
property of the profile /. 

We recall that this property is clearly true in the explicit case s = 1/2 where f(s) = (1 + 
s 2 )~( N+1 *>/ 2 . But it is not true in the limit s — > 1, i.e., in the case of the Gaussian profile of 
the Heat Equation G(s) = e~ s / 4 . Indeed, we can not obtain the same behavior for G(s) and 
sG'(s) since in this case the profile has an exponential expression. 

Proof of the proposition. We recall that 

jv 1 

(6.3) K s (x,t) = t-^f 2s {l,t~^\x\) 



(0)) where /2 S (l,x) is a continuous strictly positive function on R w of radial type, which is 
explicitly given by the expression 



sN/2 x ■ 



f 2s (l,x)= (2tt) n/2 \x\^- l / e~ u oj^ J v {\x\u)du) 



2s N 







1 



(27r)^/ 2 |x|^ Jo 



■2* 



N 



e \WJ uj^J v (uj)duj, v = {N-2)/2, 
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where J M denotes the Bessel function of first kind of order /i. For simplicity, we denote 
f(r) = /2 S (l,x), r = \x\ since /2 S (1, •) is a radial function: 

(6.4) / ( r ) = _L_ r -^^ ° e -(7) 2a a;f J v (u)du, u = (N-2)/2. 

Next, we prove an intermediate result, concerning the behavior of the derivative /'. 



Lemma 6.1 Let s 6 (0,1) and let f(r) = f2 S (l,x) be defined by (6.4). Then 

]hn r N+2s (Nf(r) + rf(r)) = - s H 2s+1 -^{^ s)T {s)T (s + ^\ 
In particular, we prove that rf'(r) ~ — r ~( N + 2s ) f r large r. 
Proof. We compute the derivative with respect to r 

Therefore 



rf{r) = -Nf(r) + r~ N J 2s (^) ' e~W '^J^du = (I) + (II) 

where (T)=—Nf(r), and (II) is given by 



(II) = Is^j, r-(^) / ( -^r^J v UUU. 



o 



According to formula (8.2), we can write 



u)Jn_-, (ui) = NJn (uj) — ojJn,, (u), 

2 2 2 

and therefore 

(II) = 2iV S -i— r-(^) /V(v)V^-V^)^ 
(2vr) A '/ 2 7 2 

1 f°° , \2s 







Then, according to Polya (see Blumenthal [7]) 



lim / e w 2s+ 2 1 Jjv (u)dw = - sin7rs / u 2s+ 2 1 KN_(uj)du. 

r->ooJ Q 2 vr Jq 2 



and 



lim / e Vr/ w s+ 2 Jjv (o;)dw = — sin7rs / oj s+ 2 X jv , -, (cj)cZu;. 

i — >oo Jq 2 ' 7T ' ' 
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Here the function are described in the paper of Erdelyi [T7] (not to be confused with 
K s (x,t)). Moreover ([T7j page 51) we have 



L 2 



N 



Therefore, 



oo 2s+ ^KN +1 (u)doj = 2 2s+ f" 1 r (s + y + 1 ) n.s). 



lim r N+2s (rf'ir) + Nf(r)) = -2sd(N, s), 



d(iV, S ) := s2 2s -J^( S in7r S )r( S )r (« + y) 



where 
(6.5) 

If we write this result as 

rN-l (rf'(r)+Nf{r)) ~ -2sCi(iV, s)r" 2s " 1 
by integrating we obtain r N f(r) ~ C\(N, s)r~ 2s , that is 

/(r)~C 1 (AT, S )r-( iV+2s ), 
which is exactly the result proved in [7]. Moreover, we obtain that 



that is 



lim r N+2s rf{r) = -(N + 2s)C 1 (N,s), 

1 — >oo 



r f'( r ) ~ ~ r ( N+2s ) for large r. 



□ 



We complete the proof of Proposition 6.1 on the behavior of the fundamental solution for 
large values of r/ = \x\ t~ l l 2s . 

Proof. The Fundamental solution is given by 

n . , 1 , ,. 
K s (x,t) =tsf(t-^\x\). 



We compute the derivative in the t variable. According to the scaling formula (6.3) we obtain 

d , \ N n -, i 1 jv_i_,, . „,. i . . 

j t K s (x,t) = --i s Vl* 5 |i|) - ^4=5 a W 

1 Y 



2s 



By Lemma 6.1 we know that 

Nf(v) + Vf'iv) ~ -2sC x (N, s)r/- (iV+2s) , for large 77, 



where Ci(N, s) is a positive constant given by formula (6.5). Therefore 



d 



#,(z,t) ~r£- 1 Ci(JV,s)7T (JV+2 ' ) = Ci(N, s)\x\~( N+2a \ for large 



□ 
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6.2 Self-similar solutions of the linear diffusion problem 

We study the existence, uniqueness and properties of self-similar solutions of the form 

(6.6) U(x,t) = t ai F(t l3l \x\) 

of the linear problem (the FPM Equation) 

{ U t + {-A) S U = for x £ M N and t > 0, 

^ ' \ U(0,x) = U {x) = C\x\^. forxGlR^, 

where C > 0, and < 7 < 2s is given. The constants ati, Pi 6R will be determined such that 



U(x,t) is a self-similar solution of Problem (6.7). 



Existence of a solution U to Problem (6.7) follows from paper |5], since the initial data 
Uo(x) = |x| 7 with 7 < 2s belongs to a suitable weighted space L 1 (R N , ipdx) . 

Let 7] = t^\x\. Then, 

U t (x,t) = ai t a ^ l F{r,) + /3 1 t a ^ 1 7 1 F'{r 1 ), 

(-A) s U(x,t) = t ai (-A) s (F(t l3l \x\)) = t a H 2 ^ s {-A) s F{r,). 
We obtain a first relation on the parameters: a\ — 1 = a\ + 2(3\s, and then j3\ = — 
Equation. The profile F satisfies the equation 

a x F{rf) + PivF'(r]) + (-A) s F(r]) = 0. 

Self-similarity condition. The equation is invariant under transformations of the form 

T x U(x,t) = \- ai U(X~ Pl x,\t). 
Therefore, u = T\u. We apply this to the initial data 

T x U(x,0) = X- ai U(X-^x,0) = X- ai -^\x\i 
and then a\ = —7 Pi- We obtain the exact value of the similarity exponents 

(0.8) „ = ft = -±. 



Notice that a± > and ft\ < 0. As a solution of the linear problem (6.7), U(x,t) can be 
computed convolution with the kernel K(-,t) 

U(x,t) = (K(-,t)*U )(x)= [ K(y,t)U (.x-y)dy. 

Since the initial data is a radial function Uq(x) = {x^ , then by the properties of the kernel K, 
U will also be a radial function, and therefore the profile F is radial. 
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Lemma 6.2 (Properties of the profile) The profile F is monotone non- decreasing and it 
satisfies r)F' < C2F, for all rj > 0. 



Proof. 

I. Monotonicity property. In order to prove the positivity of F we will make use of the 
Alexandrov Symmetry Principle and we prove that U(x,t) is radially increasing in the space 
variable x E M N . 

We start with increasing radial initial data Uq(x) = |x| 7 . We approximate Uq with a sequence 
of radially symmetric and bounded functions Uo n E L 1 (M JV ) such that Uo n (r) — > C n 7 as 
r — > 00 and vo n (r) = C rC — Uo n (r) £ L 1 (M Ar ). Let u n the solution of Problem (6.7) with 



initial datum vq u . We may apply the Alexandrov Symmetry Principle (that we explain in 
detail below) to v n to conclude that it is radially symmetric and decreasing w.r.t. the space 
variable. We then put U n (x,t) = C n 1 — v n (x,t), which is radially symmetric and increasing, 



and solves (6.7) with initial datum Uo n . We pass now to the limit n — > 00 to get the same 



conclusion for U. 

Applying the Alexandrov Symmetry Principle. We fix two points x and x' in ~M. N such 
that \x\ < \x'\. Let H denote the hyperplane perpendicular on the line xx 1 . Let Oi and VL2 be 
the two sets delimited by the hyperplane H such that the origin is contained in Oi. Let II the 
symmetry with respect to H that maps into f^. Clearly, II(x) = x', x £ fii. Then one can 
prove that for every y £ $7i \y\ < \y'\, where y' = n(y). Since VQ n is radially decreasing, we 
get that vo n (y) > vo n (H(y)), for all y G Q\. By applying the Alexandrov Symmetry Principle 



stated in Theorem 8.1 we obtain that v n (x) > v n (x'). The arguments we used can be done 



for every pair of points \x\ < \x'\, therefore v n is radially increasing. 

II. Decay at infinity. A formal computation starting from the initial data U(x,t) —> |x| 7 
as t — > gives us that ry _7 i ? (r/) — >■ 1 as rj — > 00. Therefore 

F{rf) ~ rp for large 77. 

This characterization of the profile F gives us the following spatial decay for U(x, t) for large 
times 

C 2 |x| 7 < U(x,t) < Ci|x| 7 , for large \x\t~ 1/2s . 
Moreover, we will prove the following relation between F' and F: 

\~fF(n) — 7]F'(r])\ < rp , for large r] > 0. 

As a consequence we can characterize the derivative Ut- 

U t (x,t) = r 1 " 1 ! ( 7 F(r/) - riF'(r,)) , r, = f^x]. 
Ut(x,t) ~ t _1 |a;| 7 for large values of t _1 ^ 2s |x|. 
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The first step will be to obtain a formula for the profile F{rj). Therefore 

U(x,t) = K s (x,t)*U (x)=t-& I f(t-^\x-y\)\y\^dy, z = t~&y 

JRN 

= t 2s / f(t~^x - z)\zydz. 



Since U(x,t) has the self similar form (6.6) then 

F{t~^ x )= f f{t~^ x - z )\ z \^dz = (f*U )(t~^x), VxeK w ,t>0, 
Jr n 

that is 

F(rj)= I f{r]-z)\z\^dz, Vr? > 0. 
Jr n 

Let us continue using the notations 

F(\r,\) = F( V ), f(\ v \) = f( v ). 
We fix r\ E R . Let |^| = r/ and r/ = fje for a vector e E with |e| = 1. Then 

F (v)= f(\ z \)\V ~ z V dz = V N+1 filvvDW-y^dy, z = fjy. 
Jr n Jr n 

We differentiate in f) 

F'(f)) = fj N +^ [ [(N + j)f(\fjy\) + fjyf'(\fjy\)] \e - y^dy. 
Jr n 

Therefore 

fjF'(fj)-jF(fj) = fj N+ ' y \Nf{\fjy\) +fjyf'(\fjy\)]\e-y\" f dy, z = fjy 

Jr n 

fjF'(fj)- 1 F(fj) = f)"< [ [Nf(\z\) + zf'(\z\)] e-* dz. 

We know that Nf(r) + rf'(r) Cir - ^ +2 ') for large r. Since we deal with a convolution 

we will use the information only in the sense of modulus. We fix R > such that 



Cir -(N+2s) < \ N f(j) + z fHr)\ < C2r -(JV+2,) 5 W > R 



Then 

where 
I -- 



[Nf(\z\) + zf'(\z\)] 



dz = I + II, 



\z\<R 



[Nf(\z\) + zf(\z\)] 



z 




z 




e 


1 dz<C [ 


e 


1 dz<C^ 


V 


J\z\<R 


V 





R 



,7V 
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The second term is estimated as follows (notice that R can be taken large enough to ensure 
that R> fj since the value fj was fixed at the beginning). 



H= [ [Nf(\z\) + zf'(\z\)] 

J\z\>R 



< 



\z\>R 



V 



L 



-(N+2s) ( 1 + 



< 2 / \z 

f] J J\z\>R 

2rT 7 I \z\- {N+2s ^dz, we know 7 < 2s, 
J\z\>R 



{N+2s) 

\z\>R 

(N+2s) ( JfJ 
fj 



dz 



= 2- 



1 



1 



fp 2s - 7 
We conclude that 



R? 



-2s 



I + II < C± + G 



Now, recall that r]F'(r]) — r yF{r]) = rp(I + Therefore we have proved that 
limsupfT 7 \fjF'(fj) - 7-F(r?)| < d + C 2 ^ -»• Ci, 

which means the limit is finite. Therefore, 

|7?F'(?7) — ^yF(fj) \ < fj 1 for large fj. 



□ 



7 The Reaction Problem 

As a further evidence of the influence of the tail of the data on the propagation rate, we 
consider the purely reactive problem (no diffusion) 

(7.1) u t = f(u), x£R N , t>0, 

with initial datum uq and f(u) ~ u(l — u) ~ f'(0)u. It is easy to see that when we simplify 
f{u) to f'(0)u = au, the exact solution is 

u(x,t) =u Q e f '^ t . 

Let us examine the level sets in two particular cases. 

Exponential decay. By considering initial datum of the form u${x) ~ e~ x for large \x\, 
then the solution u(x, t) satisfies a similar behavior 

u(x,t) ~ e -( x2 - at ) for large x. 

The level sets u(x, t) = constant are characterized by x = \/at + c. 
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Power decay. By considering initial datum of the form uq(x) ~ |x| ( Ar+2s ) for large then 
the solution u(x, t) is such that 

u(x,t) ~ e at \x\' (N+2s) - 

a i 

The level sets u(x,t) = constant are characterized by \x\ ~ e N + 2s . 

Conclusion: the influence of fractional diffusion: For \x\ large, the solution of the 



reaction-diffusion Problem (1.1) behaves like the solution of Problem (7.1), that is, the non- 
diffusion case. The fractional diffusion term (— A) s u does no change the basic behaviour of 
the solution for large \x\. This fact has been also observed by King and McCabe in |19| in the 
fast diffusion case with the standard Laplace operator. 

8 Appendix 



8.1 Concept of solution to Problem (1.1) 



According to |13j there exists a unique mild solution of Problem (2.1) corresponding to the 
initial datum no £ L l (K N ), < «o < 1, constructed by means of the tools of semigroup 
theory. Moreover, such n is in fact a strong solution of the equation. In the case m > 1, the 
(jot regularity of the solution follows from |5j, and this has been extended to m < 1 up to the 
extinction time (if there is one). Quantitative estimates of positivity of the solution for any 
m > corresponding to non-negative data have been proved in [8]. Recently, C 1,a regularity 
of strong solutions was proved in [14J. 

As a consequence one obtains by rather standard methods the existence, uniqueness and 



regularity properties of the solution to Problem (1.1) corresponding to the initial datum no £ 
L 1 (R N ), < n < 1. In order to prove the existence of a solution of the problem ut + L s u m = 
cu, the idea is to prove that the map no i- > v = e~ ct u is a m-w-accretive operator. Standard 
properties, like the maximum principle hold also in our setting. 

A more detailed analysis of these properties is beyond the purpose of this work. 
8.2 A version of the Maximum Principle 

We need an interesting version of Maximum Principle proved by Cabre and Roquejoffre in 
[11], Lemma 2.9, suitable for comparisons in which fractional laplacian operators are involved. 

Lemma 8.1 Let N> I, sG (0,1), < 7 < 2s. Letv G C^flO, 00); X 7 ) satisfy v(-,t) £ D^A) 
for all t > 0. Let r : [0, 00) be a continuous function and define 

Q lI = {( x ,t) e (0,oo) x R N : \x\ < r(t)}. 

Assume in addition: 
(HI) v(0,-) < inW N . 
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(H2) v<0 in ((0, oo) x R N ) \ Qj. 
(H3) a(x,t)vt + L s v < bv in Qj. 
Then v < in (0,oo) x R N . 

Although the equation we have is different, the proof as in [11] still works (with inessential 
modifications). 

8.3 Alexandrov Reflection Principle 

We recall the version of Alexandrov's symmetry principle that holds in the case of the nonlinear 
parabolic problem 

(8.1) u t = L s u m , u(0,x) = u (x), 

posed in H N , with L s = (— A) s , m > 0, uo € L 1 (M Ar ). Let us take a hyperplane H that 
divides M. N into two half-spaces Q,\ and and consider the symmetry n with respect to H 
that maps fii into ^2- The following result is proved as Theorem 15.2 in 



Theorem 8.1 Let u be the unique solution of Problem (8.1) with initial data uq. Under the 
assumption that 

uq(x) > uo(n(x)) in fix 

we have that for all t > 

u(x, t) > u(TL(x), t) for x G Oi. 
8.4 Bessel functions of first kind 

The Bessel function of first kind can be introduced through a series expansion, cf. p], 

(_l)fc fZ\ 2 *+» 



M*) = E 



^ /,! T(k + fi + l) V.2 



We mention the following recurrence formulas: 

1 

2 



J'n{z) = ~ (Jv-i( z ) ~ Jn+i{z)) , for n ^ 0. 



3-2) J^z) = —{J^ 1 (z) + J^ +1 {z)) 



4( z ) = ~ J i(z), 

z 

2/x 



(8.3) / K a {t)t h - l dt = 2 fe - 2 r ( h -±± V ), K(6±a)>0. 
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Comments and Open problems 

• There are critical values of the speed a which we do not cover in this work: o\ for m c < 
m < mi; 02 for <m\ < m < 1; respectively, (02, 03) for m > 1. The analysis of those cases 
leads to long new developments. 

• Is there a definite profile function that represents up to translation the shape of the solution 
in the region where it varies in a marked way to join the level u = 1 to the level u = 0? Maybe 
for s = 1/2 this question is easier. 

• For reasons of length and novelty, the case m < m\ is not studied. For the corresponding 
fractional fast diffusion equation there appears the phenomenon of extinction in finite time. 
King and McCabe in [19] give an idea on the asymptotics in this range of parameters. 

• A detailed numerical treatment of these problems for the case m 7^ 1 is needed, see in this 
respect [2~T] . 

• There are other interesting directions in this class of problems. Thus, in a recent paper jS], 
the authors investigate the model 

u t (x, t) + Au(x, t) = n(x)u -u 2 , x £ R N , t > 0, 

where the function \i is supposed periodic in each spatial variable Xi and satisfy < min \i < 
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